Binding of Similarly Charged Plates: A Global Analysis 



O 

o 
o 

(N 

in 

(N 



o 

CO 



o 
o 



> 

m 

OS 

o 
o 
o 

-I— > 



I 

c 
o 
o 



X 
S3 



Andr G. Moreira and Roland R. Netz 
Max Planck Institute of Colloids and Interfaces, 144^4 Potsdam, Germany 

(February 1, 2008) 

Similarly and highly charged plates in the presence of multivalent counter ions attract each other, 
leading to electrostatically bound states. Using Monte-Carlo simulations we obtain the inter-plate 
pressure in the global parameter space. The equilibrium plate separation, where the pressure changes 
from attractive to repulsive, exhibits a novel unbinding transition. A systematic and asymptotically 
exact strong-coupling field-theory yields the bound state from a competition between counter-ion 
entropy and electrostatic attraction, in agreement with simple scaling arguments. 



Experimentally, it has been known for a long time 
that highly charged planar surfaces attract each other in 
the presence of multivalent counter ions, inducing bound 
states. This electrostatic binding restricts the swelling 
of calcium clay particles |jj and leads to much reduced 
water uptake of charged lamellar membrane systems [El . 
Attractive forces between charged surfaces have also been 
observed with the surface force apparatus M. Monte- 
Carlo simulations indeed confirmed that for a given sur- 
face charge density there exists a threshold counter-ion 
valence above which attraction can be observed over some 
range of plate separations B . 

Theoretically, these observations came as a surprise, 
since the mean-field or Poisson-Boltzmann (PB) theory, 
which works usually quite well for charged systems, pre- 
dicts only repulsive forces between charged objects ||. 
This contradiction between observation and PB predic- 
tion resulted in an immense theoretical activity, which 
aimed at understanding the simple model of two uni- 
formly and similarly charged planar surfaces interacting 
across a gap of width d filled with point-like counter ions. 
Clearly, reality is much more complicated due to addi- 
tional interactions, but even this simple model, which 
we will consider in the following, is quite challenging. 
A number of approaches were proposed which incorpo- 
rate counter-ion correlations that are neglected within 
PB. The first were integral-equation theories pi, per- 
turbative expansions around the PB theory J7]|], and 
local density- functional theory B, which compare well 
with simulation results and exhibit attraction. If the 
two plates are far apart from each other, the counter- ion 
clouds can be viewed as condensed on the plates, and the 
resulting simplified model can be solved within a Gaus- 
sian [fyl or harmonic-plasmon approximation pjfl . These 
approaches either involve numerics and do not provide 
much physical insight, or they are valid for large sep- 
arations and cannot be used to characterize the bound 
state. 

Of great significance is the fact that when the elec- 
trostatic force is attractive, the equilibrium distance d* 
between the plates is smaller than the typical lateral dis- 



tance a between counter ions (as we will demonstrate 
later on), thus rendering a quasi- two-dimensional layer 
of counter ions p0[ . In the first part of this paper, we 
will use this fact and present a scaling argument for the 
attraction between two plates, valid for d <C a. We will 
then demonstrate that this scaling analysis is equivalent 
to the leading order of a systematic field-theory, valid 
in the strong-coupling (SC) limit (corresponding to large 
plate-charge density a or large counter- ion valence q), 
where it agrees with extensive Monte-Carlo (MC) sim- 
ulations. Our MC results span the complete parameter 
space. Whenever attractive forces between the plates ex- 
ist, they induce a bound equilibrium state, which exhibits 
a novel unbinding transition. 

The simple scaling argument for attraction between 
charged plates starts with partitioning the system into 
isolated counterions sandwiched between two finite plate 
segments of area A = q/2a. Neglecting ion-ion interac- 
tions should be valid for d <§; \f~A. Denoting the distance 
between the counterion and the plates as x and d — x, 
respectively, we obtain for the electrostatic interaction 
energies, in units of fcgT and for d <C yA, the results 
W\ = 2Tr£Bqc r x and W<2 — 2-K£Bqo{d — x), respectively, 
as follows from the potential at an infinite charged wall 
and omitting constant terms. The sum of the two in- 
teractions is W1+2 — W\ + W% = 2-KlBqvd which shows 
that i) no pressure is acting on the counterion since the 
forces exerted by the two plates exactly cancel and ii) 
that the counterion mediates an effective attraction be- 
tween the two plates. The interaction between the two 
plates is proportional to the total charge on one plate, 
Aa, and for d <C \f~A given by W\% = —2itA(.BG 2 d. Since 
the system is electroneutral, q = 2Aa, the total energy 
is W = W12 + W\ + W 2 = 2nA£ B cr 2 d, leading to an 
electrostatic pressure P e i = —d(W/A)/dd = ~2tt£b>J 2 . 
The two plates attract each other. The entropic pressure 
due to counter-ion confinement is P en = 1/ Ad — 2a j qd. 
The equilibrium plate separation is characterized by zero 
total pressure, Ptot — Pel + Pen — 0, leading to an 
equilibrium plate separation d* — l/n£Bq<y. By con- 
struction, the derivation for d* is valid only for d* < a 



(where a is the average lateral distance between ions 
as defined by ira 2 = q/2a), equivalent to the condition 
S = 2-K£ 2 B aq 3 > 4, i.e. for large values of the coupling 
constant S. Surprisingly, these results for P to t and d* 
become exact in the SC limit 5 — > oo, as we will demon- 
strate in the following. In fact, the attraction between 
charged plates, as derived here for d <C a, is conceptu- 
ally simpler than the PB result of repulsion, because the 
latter case involves many-body effects. 

To proceed with our systematic field theory, consider 



the partition function for N counter ions confined be- 
tween two parallel plates at distance d 
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where the Heavyside function is defined by 0{z) = 1 for 
z > and zero otherwise. Introducing the counter-ion 
density operator p(r) = Ylj=i ^( r — r j) ^ ne Hamiltonian 
can be written as 



n= e -fjdrdr'[qp(r) 



aS(z) - aS(d - z)]v(r - r')[qp(r') - aS(z') - aS{d - z')] - / dr p(r)h(r) 



(2) 



where v(r) — 1/r is the Coulomb interaction and the field h has been added to calculate density distributions 
later on. The characteristic length scales are the Bjerrum length £b — e 2 /A-neksT and the Gouy-Chapman length 
p = l/2ir£Bqcr, which measure the distance at which the interaction between two unit charges and between a counter 
ion and a charged wall reach thermal energy, respectively. Rescaling all lengths by the Gouy-Chapman length according 
to r = /if and d — pd, the Hamiltonian becomes 

H = t-^h /dfdf'[27rSp(f) - 6(z) - 6{d- z)]v(r - f')[27rEp(r') - 5{z') - 5{d- z')} - f dr p(f)h(f) (3) 

and thus only depends on the coupling parameter S = 2nq 3 £ 2 ^a. At this point we employ a Hubbard-Stratonovitch 
transformation, similar to previous implementations of a field theory for charged systems 13 1, followed by a Legendre 
transformation to the grand-canonical ensemble, Q — J^n ^ 
operator follows from Poisson's law as i; -1 (r) = — V 2 (5(r)/47r, 
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Zjy, introducing the fugacity A. The inverse Coulomb 
which leads to 



Q = 
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where we introduced the notation Z v — %/det v and 
the rescaled fugacity A is defined by A = 27rA/i 3 S = 
\/(2tt(j 2 £b)- The expectation value of the counter- 
ion density, p(r), follows by taking a functional deriva- 
tive with respect to the generating field h, p(r) = 
S\nQ/6h(r)p 3 , giving rise to 



p(r) 



P(r) 
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The normalization condition for the counter-ion distri- 
bution, p J dzp(z) = 2a /q, which follows directly from 
the definition of the grand-canonical partition function, 
leads to 



A 



di(e 



5 >) =2. 



(0) 



This is an important equation since it shows that the ex- 
pectation value of the fugacity term in Eq. (|j) is bounded 
and of the order of unity. Let us first repeat the saddle- 
point analysis, which, because of the structure of the 
action in Eq.(Q), should be valid for H <C 1. The 
saddle-point equation reads d 2 <fi(z)/dz 2 = 2iAe^ l ^ z \ 
The solution of this differential equation is i(j){z) — 

2 In cos ( A 1 / 2 ^ — d/2] j . The normalization condition 



leads to the equation A 1 / 2 tan[dA 1 / 2 /2] = 1, which is 
solved by A ~ 2/d for d < 1 and A ~ 7r 2 /4d 2 for d > 1. 
From Eq. (ra) , the rescaled density distribution of counter 
ions is given by the well-known PB result 



p(z) = 1/cos 2 U 1 / 2 [z-d/2]\. 



(7) 



Let us now consider the opposite limit, when the cou- 
pling constant S is large. In this case, the saddle-point 
approximation breaks down, since the prefactor in front 
of the action in Eq. (H) becomes small. Since the fugacity 
term is bounded, as evidenced by Eq.(H), one can expand 
the partition function (and also all expectation values) in 
powers of A/S (which is equivalent to a virial expansion) . 
For the expectation value determining the density Eq. (||) 
the leading two orders in the virial expansion are 
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The normalization condition Eq.(0) can then be solved 
by an expansion of the fugacity in inverse powers of the 
coupling constant, A = Ao + Ai/H + . . .. We obtain 
A = 2e~ v ^/ 2 jd and thus the density distribution is (in 
agreement with our scaling analysis) to leading order in- 
deed a constant given by 
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FIG. 1. MC results for the rescaled counter-ion density 
dp as a function of the rescaled distance from the wall z/d 
in the a) PB limit for S = 0.5 and in the b) SC limit for 
S = 100 for various plate separations d — d/p, = 1.5 (open 
diamonds), d = 10 (open stars), and d — 30 (open triangles). 
In a) MC results agree well with the corresponding PB pre- 
dictions (Eq.([7|), solid lines), whereas in b) results for d — 1.5 
agree with the asymptotic SC prediction, Eq.(H) (dashed line) 
and for d = 30 with a double-exponential curve (see text). 



p(z) = 2/d+0(E^) 



(8) 



In Fig. la we show counter- ion density profiles obtained 
using MC simulations |14| for small coupling parameter 
5 = 0.5 for various plate distances, which are well de- 
scribed by the PB profiles Eq.(0) shown as solid lines. 
Fig. lb shows that for 5 = 100 PB (thin solid lines) 
is inadequate |]15). As suggested by our scaling analy- 
sis, the asymptotic SC result Eq.(g) should be valid for 
d/a = d/S 1 / 2 < 1 only, since otherwise ion-ion interac- 
tions become important. For d = 3/2 (open diamonds) 
we find d/a = 0.15, and indeed Eq.(||) is accurate. For 
d = 10 (open stars) we find d/a = 1, the density pro- 
file is neither described by Eq.(B) nor (pi). Finally, for 
d = 30 (open stars) we find d/a = 3, the two layers are 
decoupled and the density profile is described by a double 
exponential p(z) = (e~ z + e z ~ d )/(l — e~ d ) (dashed-dotted 
line), which is the superposition of the density profiles of 
two isolated charged surfaces in the SC limit Jig ], The 
crossover from PB to SC is demonstrated in Fig. 2, where 
we plot density profiles for fixed plate separation d = 2 
for various coupling parameters S. 
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FIG. 2. MC results for rescaled counter-ion density pro- 
files p = p/2-k£bv 2 for fixed plate separation d — d/p = 2 
as a function of the rescaled distance z — z/p from one wall. 
Symbols correspond to coupling parameters H = 0.5 (open 
diamonds), S = 10 (filled diamonds), S = 100 (open stars), 
and E = 10 (open triangles), exhibitingclearly the crossover 
from the PB prediction (solid line, Eq.(J7])) to the SC predic- 
tion (broken line, Eq.(0)). 



Using the contact value theorem, the pressure P be- 
tween the two plates, which follows from the partition 
function via P = d\n Q\/Apdd 1 is related to the counter 
ion density at a plate by J|,[L3| 



P 



2tt£ b <7 2 



p{d) - 1. 



(9) 



The first term is the entropic pressure due to counter- 
ion confinement, the second term is due to electrostatic 
interactions between the counterions and the charged 
plates. Numerically, the contact ion density p{d) is ob- 
tained from the density profiles by extrapolation. In 
Fig. 3 we show numerical pressure data for various val- 
ues of 5. Attraction (negative pressure) is obtained for 
5 > 10. The numerical pressure for 5 = 0.5 (open di- 
amonds) agrees well with the PB prediction (solid line), 
which from Eqs.([|) and (jfy is given by P = A with A 
determined by A 1 / 2 tan^A 1 / 2 / 2 ] = 1- Tne small dis- 
tance range of most data, and the complete pressure 
data for S = 10 5 (open triangles) are well described by 
the SC prediction (broken line). It results from combin- 
ing Eqs.(W) and (||) and is given by P — 2/d — 1, from 
which the equilibrium separation, which corresponds to 
the minimum of the effective plate-plate interaction, is 
obtained as d* = 2. Incidentally, this is exactly the scal- 
ing prediction for the pressure derived in the beginning 
of this paper. 




FIG. 3. MC results for the rescaled pressure 

P = P/2tv£bo' 2 as a function of the rescaled plate separa- 
tion d — d/fi for the same parameter values as in Fig. 2 (and 
S = 20, open squares), compared with the PB prediction 
P = A (solid line) and the SC prediction P = 2/d~ 1 (broken 
line) . 

Finally, combining all pressure data, we obtain the 
global phase diagram shown in Fig. 4, featuring two re- 
gions where the inter-plate pressure is attractive and re- 
pulsive. The dividing line between those regions, which 
corresponds to the equilibrium plate separation d* in the 
bound state, is determined over 4 decades of the coupling 
constant S. In the limit of large coupling constants, the 
phase boundary saturates at d* = 2, in agreement with 
our scaling argument and the leading result of our SC the- 
ory. The threshold coupling constant to observe attrac- 
tion between charged plates is S* sa 10. As 5 — ► S* from 
above, the equilibrium plate separation diverges contin- 
uously to infinity. This constitutes a novel unbinding 
transition, which experimentally could be observed with 
charged lamellar or clay systems by raising the tempera- 
ture. 

The plate separation equals the lateral ion-distance 
on a line in Fig. 4 determined by d/a = d/S 1 / 2 = 1, 
which crosses the equilibrium-separation line at d ss 3. 
For the most part of the equilibrium-separation line in 
Fig. 4 the counter-ion distribution is therefore indeed two- 
dimensional, and many-ion effects can be neglected ex- 
cept very close to S*. Correlations between counter ions, 
except the lateral exclusion correlation which keeps ions 
apart, are therefore mostly unimportant in the bound 
state. This explains why the simple single-ion scaling ar- 
gument advanced in the beginning, which turns out to 
be exact in the limit S — > oo, works so well. More de- 
tailed comparison between MC and next-leading terms 
of our SC theory, which are equivalent to higher-order 
virial terms, together with a detailed discussion of the 
significance of Wigner crystallization (which occurs at 
5 ~ 15600 in the limit d — ► 0) will be published shortly. 
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FIG. 4. Global phase diagram featuring the regions of re- 
pulsive and attractive pressure as a function of the rescaled 
plate separation d — d/jj, and coupling strength H. The di- 
viding line denotes the equilibrium plate separation d* , which 
saturates at d* — 2 for S — > oo and which diverges as S ap- 
proaches the critical value H* ~ 10 from above. 
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